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ABSTRACT 

We have observed the giant radio galaxy Centaurus A and the globular cluster oj Centauri in the TeV 
energy region using the CANGAROO-HI stereoscopic system. The system has been in operation since 
2004 with an array of four Imaging Atmospheric Cherenkov Telescopes (lACT) with '^lOO-m spacings. 
The observations were carried out in March and April 2004. In total, approximately 10 hours data 
were obtained for each target. No statistically significant gamma-ray signal has been found above 
420 GeV over a wide angular region (a one-degree radius from the pointing center) and we derive flux 
upper limits using all of the field of view. Implications for the total energy of cosmic rays and the 
density of the cold dark matter are considered. 

Subject headings: gamma rays: search — galaxy: individual (Centaurus A) — globular cluster: indi- 
vidual {lo Centauri) 



1. INTRODUCTION 

Centaurus A (Gen A, NGC 5128, J1325-4301) is one 
of the best examples of a radio-loud AGN. Viewed at 
~60° from the jet axis (Graham et al. 1979; Dufour et 
al. 1979; Jones et al. 1996), it has been classified 
as a "misaligned" BL Lac type AGN (Morganti et al. 
1992). Estimates of the distance to Gen A range from 2 
to 8 Mpc. In this paper, we adopt a value of 3.5 Mpc (Hui 
et al. 1993). Due to its proximity and high luminosity. 
Gen A has long been considered a good TeV gamma- 
ray candidate. A detection of high energy gamma-rays 
from Gen A was reported in the 1970s, however many 
subsequent attempts have not been successful. 

First, the Stellar Interferometer, located near Narrabri, 
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reported a positive detection with a fiux of /(> 
0.3 TeV) ^ (4.4±1) X 10"" photon cm-^sec-i (Grindlay 
et al. 1975). 

At higher energies, the Buckland Park array and the 
JANZOS Observatory also reported the detection of 
gamma-rays. The Buckland Park flux, measured be- 
tween 1984 and 1989, was /(> lOOTeV) ~ (7.4 ± 2.6) x 
10" 1^ photon (Clay et al. 1994) and the JANZOS flux, 
during the period Apr-Jun 1990, was /(> llOTeV) ^ 
(5.5 ± 1.5) X 10^12 piiotoncm^^sec-i (Allen et al 1993). 
Given the signiflcant attenuation expected in the gamma- 
ray flux during interactions with the cosmic microwave 
background at these energies, these measured fluxes im- 
ply much greater intrinsic fluxes (Protheroe 1986). 

However, CANGAROO-I (RoweU et al. 1999), JAN- 
ZOS (long-term) (Allen et al 1993,b), and Durham (Gar- 
ramiiiana et al. 1990) observed the Gen A nuclear re- 
gion and set upper limits on the emission in the VHE 
range. Recently, the H.E.S.S. group observed Gen A for 
4.2 hours in 2004 (Aharonian et al. 2005c). Their up- 
per limit was 5.68 x 10"^^ photon cm~^ sec~^ (1.9% of 
the Grab flux) at 190 GeV. They only investigated the 
region close to the center of Gen A, i.e., the inner jet re- 
gion. Gen A has a large structure revealed by the radio 
observations, with inner lobes extending over '^10 ar- 
cmin, a middle lobe extending over ~1 degree, and outer 
lobes extending over more than 5 degrees (Burns, Feigel- 
son, & Schreier 1983). The fluxes obtained in those 
observations are plotted in Fig. 1. The details of the 
observations are listed in Table 1. 

Gen A has displayed pronounced variability at X-ray 
energies over the last 35 years, with the Grindlay et al. 
detection coinciding with the peak of the X-ray flux over 
this period (see, e.g.. Turner et al. 1997). Data from 
the ROSAT An Sky Monitor (ASM 2007), indicate that 
in recent years, including the period of the H.E.S.S. ob- 
servations. Gen A has been near its historical minimum 
X-ray state (MJD from 53111 to 53113). 
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TABLE 1 
Details op past observations. 



mark^ 


reference 


epoch 


threshold 


significance 


method 


G 


Grindlay et al. (1975) 


1972-1974 


300 GeV 


4.5 <T 


Cherenkov'' 


J3 


Allen et al (1993) 


Apr-Jun 1990 


110 TeV 


2%d 


Shower Array 


J2 


Allen et al (1993) 


1987-1992 


110 TeV 


95% CL. UL. 


Shower Array 


B 


Clay et al. (1994) 


1984-1989 


100 TeV 


99.4% CL. 


Shower Array 


Jl 


Allen et al (1993b) 


1988 & 1989 


1 TeV 


95% CL. UL. 


Chaxenkov'' 


D 


Carraminana et al. (1990) 


1987 & 1988 


300 GeV 


3 a UL. 


Cherenkov'' 


CI 


Rowell et al. (1999) 


Mar-Apr 1995 


1.5 TeV 


3 cr UL. 


lACT 


H 


Aharonian et al. (2005c) 


Apr 2004 


190 GeV 


99% CL. UL. stereoscopic lACT 



''Mark in Fig. 1. 

'^No imaging Cherenkov telescope. 
"^Spectral index consistent with 7=-1.7. 
'^Chance probability. 




Energy [TeV] 

Fig. 1. — Summary of the past observations. The filled marks 
are the positive detection and the open upper limits. The details 
of these observations were listed in Table 1. The dashed line is 
a spectrum proportional to with H.E.S.S.'s upper-limit level 
flux. 



CANGAROO-III has data outside of this minimum pe- 
riod (MJD from 530080 to 530107). We also report the 

results over a wide region. Part of the southern outer 
lobe of Ccn A was located within the field of observations 
made of w Centauri [lo Cen). We, therefore, include the 
results of the latter observation in this report. Located 
at a distance of 4.9 kpc, oj Cen is one of the oldest, heav- 
iest globular clusters in our Galaxy. As globular clusters 
have been found to contain large numbers of millisec- 
ond pulsars, they arc interesting objects for high-energy 
gamma ray observations. 

2. CANGAROO-III STEREOSCOPIC SYSTEM 

The use of imaging atmospheric Cherenkov telescopes 
(lACTs) was established with the statistically unassail- 
able detection of the Crab nebula at TeV energies by 
the Whipple group (Weekes et al. 1989). This technique 
enables TeV gamma-rays to be selected from the huge 
background of cosmic rays with the use of the "image mo- 
ments" of the shower images (Hillas 1985). Stereoscopic 



observations, which allow the signal-to-noise ratio to be 
significantly improved, were pioneered by the HEGRA 
group (Aharonian et al. 1999). The H.E.S.S. group 
has recently reported the detection of faint gamma-ray 
sources with an angular resolution as fine as a few arc- 
minutes (Aharonian et al. 2005). 

CANGAROO-III is one of two major lACTs located 
in the southern hemisphere. The CANGAROO-III 
stereoscopic system consists of four imaging atmospheric 
Cherenkov telescopes located near Woomera, South Aus- 
tralia (31°S, 137°E). Each telescope has a 10 m diameter 
segmented reflector, consisting of 114 spherical mirrors 
made of FRP (Kawachi et al. 2001), each of 80cm diam- 
eter, mounted on a parabolic frame with a focal length 
of 8m. The total light collecting area is 57.3 m^. The 
first telescope, Tl, which was the CANGAROO-II tele- 
scope, is not presently in use due to its smaller field of 
view and higher energy threshold. The second, third, 
and fourth telescopes (T2. T3; and T4) were used for 
the observations described here. The camera systems for 
T2, T3, and T4 are identical and their details are given 
in Kabuki et al. (2003). The telescopes are located at 
the east (Tl), west (T2), south (T3) and north (T4) cor- 
ners of a diamond with sides of ~100m (Enomoto et al. 
2002b). 

3. OBSERVATIONS 

The observations were carried out in the period from 

2004 March 16 to April 19 using "wobble mode" in which 
the pointing position of each telescope was shifted in dec- 
lination between ±0.5 degree every 20 minutes (Daum 
et al. 1997) from each target: (RA, dec [J2000]) = 
(201.365°, -43.019°) for Cen A, (RA, dec [J2000]) = 
(201.691°, -47.477°) for uj Cen. One of the reason why 
we took "wobble" observation is to enlarge the effective 
FOV, the other is to average the responses of individual 
pixels. We, therefore, took LONG OFF source run of 
"wobble" mode for background subtractions in the later 
analysis. 

Data with GPS time stamps were recorded for T2, T3 
and T4 individually when more than four photomultiplier 
(PMT) signals exceeded 7.6 photoelectrons (p.e.). In the 
offline analysis stage we combine all these data when the 
three telescope's GPS times coincide. The typical trigger 
rate was 11 Hz for three-fold coincidences. Each night 
was divided into two or three periods, i.e., ON-OFF, 
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OFF-ON-OFF, or OFF-ON observations. ON-source 
observations were timed to contain the meridian passage 
of the target. On average the OFF source regions were 
located with an offset in RA of +30° or -30° from the 
target. The OFF-source observations were also made in 
wobble mode. One day was dedicated to ON and OFF 
source observations of one target, with the following day 
dedicated to the other target. The images in all three 
telescopes were required to have clusters of at least five 
adjacent pixels exceeding a 5p.e. threshold (three- fold 
coincidence). The event rate was reduced to ~7.5Hz by 
this criterion. Looking at the time dependence of these 
rates, we can remove data taken in cloudy conditions. 
This procedure is the same as the "cloud cut" used in 
the CANGAROO-II analysis (Enomoto et al. 2002). The 
effective observation times for ON and OFF source ob- 
servations were 639.4 and 586.9 min for Cen A, and 600.8 
and 429.4 min for oj Cen. The mean zenith angles were 
17.4° and 20.6°, respectively 

The light collecting efficiencies, including the reflectiv- 
ity of the segmented mirrors, the light guides, and the 
quantum efficiencies of photomultiplicr tubes were moni- 
tored by a muon-ring analysis (Enomoto et al. 2006) with 
the individual trigger data in the same periods. The light 
yield per unit arc-length is approximately proportional 
to the light collecting efficiencies. The average ratios 
of these at the observation period with respect to the 
mirror production times (indicating the amount of dete- 
rioration) were estimated to be 70, 70, and 80% for T2, 
T3, and T4, respectively. The measurement errors are 
considered to be at less than the 5% level. Deterioration 
is mostly due to dirt and dust settling on the mirrors. 

4. ANALYSIS 

Most of the analysis procedures used are identical with 
those described in Enomoto et al. (2006b). As a full 
instrumental description was given in Enomoto et al. 
(2006), we omit a detailed discussion here. There are 
some improvements in the analysis procedure from that 
in the previous paper and so we concentrate on those 
points here. 

At first, the image moments (Hillas 1985) were cal- 
culated for the three telescopes' images. The gamma- 
ray's incident directions were determined by minimizing 
the sum of squared widths (xo^ weighted by the photon 
yield) of the three images seen from the assumed position 
(fitting parameter) with a constraint on the distances 
from the intersection point to each image center (Dip). 
The Dip can be estimated from the ratio Length/Width. 
The prediction curve (/) and its error (a) were estimated 
using Monte-Carlo simulations. The constraint is; 

2 _ [Dip ~ f (Length/Width)]^ 
Xc - ^ • 

In order to balance dimensions between Xa ^-nd x^, we 
need to multiply n- < WidtK^ > by Xc: where n has 
the units of number of photo-electrons (p.e.) and is op- 
timized by Monte-Carlo simulations (n=10 p.e. for this 
analysis). Finally the minimizing variable is; 

= Xo + < Width^ > -xl 

For large images, the former term dominates and for 
small images (for example, events at large zenith angles), 



the latter dominates. With this fit, we can improve uni- 
formity of 6^ resolution especially with respect to zenith 
angle. 

In order to derive the gamma-ray likeliness, we used 
the Fisher Discriminant (hereafter FD) (Fisher 1936; 
Enomoto et al. 2006). Input parameters were 

P = (W2, W3, WA, L2, Li, LA), 

where W2,W3,WA, L2, L3, LA are energy corrected 
Widths and Lengths for the T2, T3, and T4. FD has a 
small dependence on the zenith angle, i.e., as the zenith 
becomes larger, images become smaller {FD larger). We 
corrected for this using Monte-Carlo simulations. 

We rejected events with any hits in the outermost layer 
of the cameras ("edge cut"). These rejected events cause 
finite deformations especially in the Length distribution 
which results in deformations of the FD. In this anal- 
ysis, we allowed less energetic hit pixels on the "edge" 
layer if their pulse heights were less than that of highest 
15 pixels. This cut was improved in order to increase 
acceptance. 

Since we have FD distributions for OFF-source data 

and the Monte-Carlo gamma-ray events, we can as- 
sume these are background and signal events respec- 
tively. Note that in the gamma-ray simulations we used 
a spectrum proportional to E'^ where 7=— 2.1. We can 
therefore fit the FD distribution of ON with the above 
emulated signal and real background functions, to de- 
rive the number of signal events passing the selection 
criteria. With this fit, we can determine the gamma-ray 
excess without any positional subtractions. This is a one- 
parameter fitting with the constraint that sum of signal 
and background events corresponds to the total number 
of events. These coefficients can be derived exactly ana- 
lytically. 

5. RESULTS 

The morphology (significance map) is shown in Fig. 2. 

The FOV was binned to 0.2 degree x 0.2 degree 
squares. The FD distribution inside each bin was made. 
The signal function (for gamma rays) was made using 
the Monte-Carlo simulation. The background (for pro- 
tons) function was made collecting OFF events inside 
0.6 X 0.6 degree^ centered around the corresponding sig- 
nal bin. The one-parameter x^ ht was carried out and 
the excess counts were obtained inside a one-degree circle 
(the dot-dashed circle) centered at the average pointing 
position. The pointing centers are shown by the crosses. 
At one degree radius, the acceptance decreases by 30%. 
There is, therefore, not any statistically significant ex- 
cess anywhere. The positive and negative fluctuations 
roughly agree in the Cen A field of view. In that of 
u) Cen a positive offset was observed (~0.4 a), however, 
it is not statistically significant. The peak located at 
north of cj-Ccn field and its significance is 1.5 cr. The 
point spread hmction (PSF) is shown by the dashed circle 
at left-lower place. The solid contours are obtained via 
Skyview (2007) and are the GB6 radio data (4850 MHz) 
showing the inner lobes and middle lobe, and the dotted 
contours are 408 MHz radio data, roughly showing the 
outer lobes of Cen A from Fig. 11 of Burns, Feigelson, & 
Schreier (1983). 
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Fig. 2. — Significance map. The angular bins are 0.2 X 0.2 
degree^ squares. This was obtained by the fitting procedure de- 
scribed in the text. For the background function of the fit, the 
OFF region of a 0.6 X 0.6 deg^ square centered on (but not includ- 
ing) the corresponding bin was used. The smoothing was carried 
out averaging the neighboring bins. The crosses are the pointing 
centers. Invents within one degree radius circles are plotted (the 
dot-dashed circle). The point spread function (PSF) is shown by 
the dashed circle. The solid contours are 4850 MHz radio data 
of the inner lobes and middle lobe and the dotted contours are 
408 MHz observations of the outer lobes. 



Fig. 3. — Definition of search regions. Region 1 contains the jet 
and inner lobes of Cen A and region 2 is Cen. These sizes are 
identical to the PSF. Region 3 is the middle lobe of Cen A, centered 
at (RA, dec [J2000]) = (201°.95, -42°. 85). The size is identical to 
the PSF. Region 4 is within a one-degree circle from the center of 
Cen A not including regions 1 and 3: which contains portions of the 
outer lobes of Cen A close to the nucleus. Region 5 is the overlap 
of two one-degree circles centered at (RA, dec)=(205°, —47°) and 
uj Cen: the southern extremity of the outer lobe of Cen A (Burns, 
Feigelson, & Schreier 1983). Region 6 is within one-degree circle 
from the center of ui Cen not including regions 2 and 5. 



In order to derive a UL for a point source, we choose 
0^=0.06 degree^ from the center of Cen A for the cut 
position. The 2a UL was obtained to be 39.8 events. 
The energy threshold for this analysis was obtained to 
be 424 GeV by the Monte-Carlo simulation. For w-Cen, 
a 2ct UL of 32.3 events was obtained. 

In order to derive flux ULs for the jet and lobes of 
Cen A, and lo Cen, we defined the search regions as 
shown in Fig. 3. Region 1 contains the jet and inner 
lobes of Cen A, region 3 contains the middle lobe, and 
regions 4 and 5 contain portions of the outer lobes, while 
region 2 contains w Cen. There is no high-energy astro- 
nomical counterpart in region 6 (although low level ra- 
dio emission from the southern outer lobe of Cen A may 
extend into this region). The flux upper limits for vari- 
ous regions were obtained from the FD distributions of 
the ON-source runs for corresponding areas. The OFF- 
source runs were used for the background FD functions 
and the Monte-Carlo simulations for gamma rays for the 
signal. We do not observe any statistically significant 
excesses in any regions. 

Thus we can obtain fiux ULs for region 1. Those of 
Cen A center region are summarized in Table 2. ULs are 

TABLE 2 

Flux Upper Limits above various energy thresholds for 
THE Cen a ,iet 



Excess Upper Limit Energy Threshold Flux Upper Limit 
Events GcV cm s 

308 424 0.491 X 10"" 

12.6 2074 0.903 x IQ-^^ 

5.0 6202 0.485 x 10"" 



presented for three different energy thresholds. 

Exactly the same procedure was carried out for the 
u) Cen region. Again there is no statistically significant 
excesses. The flux ULs for region 2 are summarized in 
Table 3. The flux ULs for the regions 3, 4, 5, and 6 are 

TABLE 3 

Flux Upper Limits above various energy thresholds for 
uj Cen 



Excess Upper Limit Energy Threshold Flux Upper Limit 


Events 


GeV 


^ -1 
cm s 


32.3 


471 


0.355 x 10"" 


7.9 


2216 


0.519 X 10^12 


7.0 


6518 


0.602 X 10-" 



summarized in Table 4. Here, the ULs are divided by the 
TABLE 4 

Flux Upper Limits for the Other Regions 

Region Excess Upper Limit Energy Threshold Flux Upper Limit 
Events GeV cm s ^Sr 



3 37.3 424 0.875 x 10"^ 

4 132.6 424 0.193 x lO"'' 

5 203.7 471 0.106 x IQ-'^ 

6 304.3 471 0.436 x 10"" 



solid angles of the observations, because the inner lobes, 
middle lobe, and outer lobes are candidates for diffuse 
emission. 

6. DISCUSSION 
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First, our flux ULs together with that of H.E.S.S. 
are quite low compared to the past indications of TeV 
gamma rays and also to the past UL measurements, 
demonstrating the strength of the stereo-IACT tech- 
nique. 

Next we discuss the total cosmic-ray energy and the 
cold dark matter (CDM) density for each correspond- 
ing astronomical object. Here, we assumed distances of 
Ccn A and lo Ccn to be 3.5 Mpc and 6kpc, respectively. 
For the total energy of the cosmic rays, we asstime that 
the origin of the TeV gamma-rays is the inverse Compton 
(IC) scattering of cosmic-ray electrons with the cosmic 
microwave background (CMB). The electron/proton ra- 
tio of our Galaxy is thought to be 0.1~1%. Therefore, 
10^^ erg for the electron cosmic-ray component might be 
a standard value for typical galaxies. Here, we assumed 



the electron energy spectrum of oc -E 



We in- 



tegrate the electron spectrum greater than 1 GeV in or- 
der to calculate the total energy of electron component. 

The interstellar matter densities arc not well measured 
in the corresponding regions. The electron assumption is 
reasonable, because the CMB density is well defined and 
can be considered to a lower limit for ambient photon 
density. 

Gen A is considered to be a low-energy peaked BL Lac 
object (LBL) rather than high-energy peaked object 
(HBL). The spectral energy distribution (SED) is plot- 
ted in Fig. 4. Even if a cosmic ray flux similar to our 
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Fig. 4. — SED for the jet region of Cen A. The open circle is 
the prediction of Bai & Lee, (2001). The blaclc square with the 
arrow is the H.E.S.S. upper limit (Aharonian ct al. 2005c). The 
open diamonds arc our upper limits. The expected gamma-ray 
yield with a total electron energy of 10^** erg with a 100 TeV cutoff 
energy is shown by the dashed curve. The details of the marked 
observations were listed in Table 1. 

Galaxy's exists, there is no contradiction with any past 
measurements in the all wave length region (NED 2007). 
This, therefore, is not denying higher energy component 
such as 100 TeV even in the higher level compared to our 
Galaxy, for example 10^^ erg of electrons. Actually, Bai 
& Lee, (2001) predicted a time-dependent huge energy 
flow in the TeV energy region which is shown by the open 
circle in Fig. 4. The H.E.S.S. UL is shown by the open 



square with arrow and our ULs by the open diamonds. 
These are orders of magnitude lower than the predicted 
value. However, as our measurement periods correspond 
to a quiet phase of Cen A (from ASM data), there is no 
direct contradiction of this idea. 

We rather concentrate ourselves to the quiet, stable, 
and average states of Cen A. If a high energy compo- 
nent of electrons oc E~ 
of the total energy of 10 
times stronger than that of our galaxy, the expected TcV 
gamma-ray flux is the dashed curve in Fig. 4. That is 
higher than our ULs, i.c, wc can derive a meaningful UL 
to the total electron densities. The 2a ULs (the dot- 
ted curve) for the above assumption were derived and 
shown in Fig. 5. The vertical scale is the logarithm of 



e looTev exists at the level 
'^^ erg which is one-hundred 
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Fig. 5. — One- and two-cr upper limits to the total energy of the 
electron cosmic ray flux inside the jet region of Cen A versus the 
maximum accelerated energy. The one-u line is the solid curve and 
the two-(T line is dotted. 

total electronic energy normalized to 10^° erg and the 
horizontal scale is the maximum accelerated energy of 
the electron in TcV. At 100 TcV, the limit is SxlO'^^'^erg 
which is thirty times bigger than that of our Galaxy. Al- 
though Cen A is highly active especially at radio, IR, and 
MeV energies compared to our Galaxy, any high energy 
cosmic-ray components are not detected; Cen A is a mys- 
terious object. If we assume that the target photons are 
of those IR temperature blackbody radiation, the limits 
of the total energy will be signiflcantly reduced. The fu- 
ture Gherenkov telescope array experiment (GTA) (CTA 
2007) is promising to help resolve this, as it aims to be 
two orders of magnitude more sensitive than the present 
lAGTs. If GTA does not detect any signal from this 
region, we might need to totally revise the idea of the 
cosmic-ray origin. 

According to the ASM data (ASM 2007), the X-ray ac- 
tivity of Cen A was at a minimum in the M.ID range be- 
tween 53111 53113. The H.E.S.S. observation was com- 
pletely coincident with this period, and ours partially. 
We, therefore, can divide the observational period into 
two, i.e., coincident with the ASM minimum, and the 
remaining data. The ULs for these two periods are sum- 
marized in Table 5. There is however no excess in the 
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TABLE 5 

Flux Upper Limits in two epochs for Cen A jet 



Excess Limit Energy Threshold Flux Limit MJD-53000 
Events GeV cm s day 

57^6 424 0.128 x 10"^" 80, 82, 89, 91, 107 

26.7 424 0.737 x 10"^^ 111, 113 



period of the (slightly) higher state. 

Lu Cen is the biggest and oldest globular cluster inside 
our Galaxy. There is no ionizing gas which is considered 
to be an origin of cosmic rays. We, therefore, do not 
expect any TeV radiation from this object. However, to 
know the minimum level of cosmic-ray density is impor- 
tant. There may be a lot of milli-second pulsars as has 
been found in other globular clusters. The object may 
still hold a CDM halo (Peebles 1984), as these are not 
clearly denied at the present level of measurements. The 
ULs for the total energy of electron cosmic ray compo- 
nent is derived in the same way as described previously. 
Under the assumption of a maximum electron acceler- 
ated energy of 100 TeV, the total electron energy is less 
than 4 x 10^'' erg, i.e., less than one-supernova level, 
a surprisingly quiet object considering its total mass of 
- 5 X IOHIq (van de Ven et al 2006; Meylan et al. 1995; 
Richer et al. 1991). 

Cen A is a massive object and uj Cen a much higher 
density object than our Galaxy as a whole. The upper 
limit study of the CDM density would provide meaning- 
ful information as to the mechanism of their formations. 
However note that for the case of globular clusters, CDM 
was considered to be stripped off via tidal force in the pe- 
riodic motions inside the parasite galaxy. Peebles (1984) 
discussed the existence of the dark matter in this kind of 
object. The CDM search was carried out following the 
procedure of Enomoto et al (2003). Here we additionally 
used the EGRET UL at 189 MeV (NED 2007). This is 
effective around the sub- TeV CDM particle mass region. 
Assuming the CDM particle mass of 1 TeV, the 2-a up- 
per limit for the CDM density is 1.2 Mq/pc^. These are 
quite high compared to our local density estimation of 
CDM of O.OlMo/pc^. The mass of Cen A was estimated 
to be several x 1O"M0 (Mathieu et al. 1996; Peng et al. 
2004a,b; Karachentsev et al. 2007; Woodley 2006). In 
order to derive the total mass of the CDM, we multiply 
by a factor of 1.42 x 10^'' pc'' considering our angular 
resolution of 6*^=0.06 degree^ with an assumption of the 
distance of 3.5 Mpc (Hui et al. 1993). Then the total 
CDM mass inside the volume is less than 2 x 10^^ Mq, 
i.e., a meaningless result. 

For UJ Cen, we can carry out the same discussion and 
the results are shown in Fig. 6. The result is one order of 
magnitude higher than that of Cen A. This is due to dis- 
tance of 6kpc compared to that of Cen A of 3.5 Mpc 
while the volumes are determined by the same angu- 
lar resolution of the observation. However, multiplying 
the volume factor of 4x10^ pc^, then the most CDM- 
particle-mass region is under its total gravitational mass 
of 5x10^ Mq. In most regions of CDM particle mass 
between 100 GeV and 10 TeV, we can reject the hypoth- 
esis that the gravitational mass is dominantly occupied 
by the CDM particles, they might be absorbed by our 
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Galaxy. Presently, high energy physicists believe that 
the; CDM particle mass lies inside the above region. The 
future CTA project (CTA 2007) would reveal the fact 
in a wider mass region, or the near-future LHC project 
might give a conclusive result on this subject. In cos- 
mological life, the larger eats the smaller but its core 
remains. 

7. CONCLUSION 

In this paper, we present the results of stereoscopic 
observations of Centaurus A and w Centauri with the 
CANGAROO-III telescopes. The observation period was 
2004 March 16 to April 19 and the total observation times 
were 640 min for Cen A and 600 min for w Cen. The ob- 
servation of Cen A was carried out over a similar period 
to that of H.E.S.S. The w Cen observations are the first 
reported trial by Imaging Atmospheric Cherenkov Tele- 
scopes. The analysis was carried out inside a one-degree 
(radius) circle from the average pointing position. We 
derived flux upper limits for regions containing the jet 
and inner lobes, the middle lobe, and portions of the 
outer lobes of Cen A, and center of lo Cen. The Cen A 
upper limits are, as with the H.E.S.S. limit, an order of 
magnitude lower than previous measurements. 

The upper limits for the total energy of the electron 
components of cosmic rays were calculated under the 
assumption of Inverse Compton Scatterings on the cos- 
mic microwave background, with limits of 3x10^"^ and 
4x10^'' erg obtained for Cen A and uj Cen, respectively. 

Finally we gave upper limits to the density of Cold 
Dark Matter (CDM). Around the TeV region, we ob- 
tained upper limits of its density of 2 nQpc^"* for Cen A 
and 100 Mqpc^"^ for to Cen. The limit for Cen A was 
greater than its gravitational mass, however, that for 
LO Cen was less than it. 
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